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Abstract. We present detailed spatially-resolved spec- 
troscopy results of the observation of Abell 1835 using 
the European Photon Imaging Cameras (EPIC) and the 
Reflection Grating Spectrometers (RGS) on the XMM- 
Newton observatory. Abell 1835 is a luminous (10 46 ergs 
s -1 ), medium redshift ( z — 0.2523), X-ray emitting cluster 
of galaxies. The observations support the interpretation 
that large amounts of cool gas are present in a multi-phase 
medium surrounded by a hot (kT e =S. 2 keV) outer enve- 
lope. We detect 0 VIII Lya and two Fe XXIV complexes 
in the RGS spectrum. The emission measure of the cool 
gas below kT e — 2.7 keV is much lower than expected from 
standard cooling-flow models, suggesting either a more 
complicated cooling process than simple isobaric radiative 
cooling or differential cold absorption of the cooler gas. 

Key words: Galaxy Clusters: Individual: Abell 1835 - 
Galaxies: Cooling Flows - Techniques: Spectroscopic 


1. Introduction 

Previous soft X-ray imaging and medium spectral resolu- 
tion studies have implied the presence of cool gas near the 
centers of clusters (See e.g. Sarazin 1988, Fabian 1994). 
For the Centaurus, Perseus, Virgo, and Coma clusters 
there is direct spectroscopic evidence for cool gas from the 
Einstein Focal Plane Crystal Spectrometer (Canizares et 
al. 1979, 1982, Mushotzky & Szymkowiak 1988). Based 
on simple, robust estimates for the electron densities and 
temperatures, it is estimated that in the large majority of 
clusters, gas in the central volume should cool radiatively 
to very low temperatures on a timescale short compared to 
the characteristic age of the cluster. Due to the consequent 



Fig. 1. Smoothed logarithmic intensity EPIC-MOS image of 
Abell 1835. The image is approximately 23 by 13 arcminutes. 


lack of central pressure support, gas has to continuously 
flow inward into the cluster potential. 

This picture of cluster cooling flows is incomplete in 
that the fate of the cooling gas is unknown. Material that 
has cooled to below X-ray emitting temperatures generally 
fails to be detected by emission at any other wavelengths, 
although more recently there are indications that cold gas 
may be distributed throughout the cooling flow. This gas 
betrays its presence through X-ray absorption of the cool- 
ing flow emission (but see below for evidence that such 
absorption is not confirmed at higher spectral resolution). 
A simple single phase model for the cooling gas based on 
the conservation laws for the global flow predicts a sharply 
peaked density distribution that is not consistent with the 
observed distributions, and this fact has given rise to the 
idea of a diffuse, multi-phase cooling flow (Johnstone et 
al.1992). 

Abell 1835 is a prototypical strong cooling-flow duster 
with an inferred mass deposition rate of 1750 ± 550 (Allen 
et al. 1996). The central cooling time (the time it would 
take the gas to radiate away all of its kinetic energy) is 
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less than 1 Gyr, whereas the cluster could remain rela- 
tively undisturbed by mergers for ~ 10 Gyr. So roughly 
the inner 350 kpc of the cluster should have had sufficient 
time to cool below X-ray temperatures. Allen et al. (1996) 
found Abell 1835 to be one of the strongest cooling-flows 
from the deprojection of ROSAT HRI images. The ASCA 
spectrum was consistent with a kT e = 8.5 keV collisionally 
ionized plasma, but cooling-flow spectral models were un- 
constrained (Allen et al. 1996). Allen and Fabian (1998) 
also carried out a joint ROSAT-ASCA study and found 
evidence for a massive cooling flow, as well as a relatively 
high metallicity. Until this observation, cooling-flow prop- 
erties have largely been inferred by relatively short cooling 
times. Here we present spectroscopic results on the prop- 
erties of the cool gas. 

As part of the XMM-Newton (Jansen et al 2000) per- 
formance verification program, Abell 1835 was observed 
for 60 ks. Data were obtained simultaneously with the two 
Reflection Grating Spectrometers (RGS) and the three 
European Photon Imaging Cameras (EPIC). The RGS is 
capable of achieving high spectral resolution for the entire 
soft wavelength band (5 to 38 A) for moderately extended 
X-ray sources. The spectral resolution for sources of char- 
acteristic angular size, 0, is roughly, AA « 0.1 A (6/1 1 ). 
The high throughput of the XMM mirrors allows the 
EPIC-pn and EPIC-MOS detectors to perform detailed 
spectral-spatial studies of clusters since the spatial extent 
of most clusters is resolved. Both sets of instruments can 
be used together to produce a spectral-spatial model of 
the X-ray emission. 

Throughout this paper we assume Hq = 
70 km s~\ Mpc -1 , D m = 0.3, and Oa — 0.7. This 
implies an angular distance of 1280 Mpc and a luminosity 
distance of 1600 Mpc for Abell 1835. One arcsecond 
corresponds to 6.2 kpc. Errors are quoted at the 90% 
confidence level. 

2. Overall Temperature Structure 

In order to study the cooling flow properties of Abell 1835 
in detail, we first use the spatially-resolved EPIC-MOS 
and EPIC-pn spectra to characterize the temperature and 
density structure. The EPIC instruments and their cali- 
bration are described in Turner et al. (2000) and Struder 
et al. (2000). An image of Abell 1835 is shown in Fig. 1. 

We performed spectral fits with EPIC-MOS and EPIC- 
pn on two concentric annuli in Abell 1835. EPIC-pn fits 
axe shown in Figs. 2 and 3. Background spectra were used 
from Lockman Hole observations (Revolution 70 and 73, 
lOOks observation) or off-source regions and the spectra 
were corrected for energy-dependent telescope vignetting. 
(See also Arnaud 2000) The spectral and spatial responses 
do not vary strongly over the chosen annuli. Periods with 
high background rates were removed. We fit the spectra 
to single temperature models with the XSPEC vll pack- 
age (Arnaud 1996) using the MEKAL (Mewe, Kaastra, & 


Table 1. Basic Temperature Structure: EPIC-MOS Spectral 
Fits 


Parameter 

e < o.5' 

0.5' < e < 2' 

kT (keV) 

6.0 ±0.3 

8.2 ±0.4 

^Fe 

0.35 ± 0.05 

0.28 ±0.05 

7V H (10 2O cm -2 ) 

2.6 ±0.4 

2.1 ±0.4 

* 2 /dof 

1.08 

0.97 


Table 2. Basic Temperature Structure: EPIC-pn Spectral Fits 


Parameter 

6 < 0.5' 

0.5' < e < 2' 

kT (keV) 

5.7 ±0.2 

8.4 ±0.4 

-dp e 

0.34 ± 0.05 

0.28 ± 0.06 

Ar H (10 2O cnr 2 ) 

1.8 ±0.3 

1.2 ±0.3 

X 2 /dof 

1.11 

0.98 


A 1835 M0S1 +MOS2 + pn R<0.5’ 



Energy (keV) 

Fig. 2. EPIC-pn (green), EPIC-MOS 1 (black) and EPIC- 
MOS 2 (red) spectral fit and residuals of the inner region of 
Abell 1835. 

Liedahl 1995) spectral model. Elemental abundance ratios 
were held fixed to their solar values when they were not 
well-constrained. 

Within a 0.5 arcminute radius circle, Abell 1835 has 
a temperature of ~ 6 keV and an iron abundance of 0.35 
relative to solar. The outer region, 0.5 < r < 2 arcminutes, 
has a temperature of ~ 8 keV and an iron abundance of 
0.28 of solar values. Column densities were consistent with 
the galactic value in the MOS fits and slightly lower than 
galactic in the PN fits. There is no evidence for a strong 
abundance gradient. The results of the fits are summarized 
in Tables 1 and 2. A combined fit gave a silicon abundance 
of 0.19 ± 0.20 and oxygen abundance 0.33 ± 0.15. 

The spectral fits are consistent with the interpretation 
that there is a large volume in the center of Abell 1835 
with cooler gas. We use the single temperature fits to the 
outer annulus to provide an upper temperature limit for 
the RGS analysis. 
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Fig. 3. EPIC-pn (green), EPIC-MOS 1 (black) and EPIC- 
MOS 2 (red) spectral fit and residuals of the outer region of 
Abell 1835. 

3. Cooling-Flow Spectroscopy 

By studying the RGS spectrum, we can measure both the 
relative emission measure of cold gas and study the spec- 
tral properties of the cooling flow in detail. 

The properties of the RGS and its calibration are de- 
scribed by den Herder et al. (2000). The RGS instrument 
effectively works as a long-slit spectrometer, providing 
spectral information in one direction (somewhat degraded 
by the spatial extent of the source) while providing 15 
arcsecond imaging information in the cross-dispersion di- 
rection. For each photon, the focal plane detectors record 
the position along the dispersion direction (which gives the 
dispersion angle), the position in the cross-dispersion di- 
rection, and the CCD pulse-height. An energy dependent 
joint dispersion coordinate/pulse-height filter separates 
the spectral orders. In this analysis, we jointly use first 
and second orders. For all spectra presented here, we as- 
sign wavelengths based on the nominal position of the cen- 
ter of the X-ray emission (a = 14 /l l m l # , S = +2°52'43'). 

The spectrum for both RGS instruments is shown in 
Fig. 4. The spectrum has been corrected for redshift, di- 
vided by the effective area and exposure, and a back- 
ground model is subtracted for display purposes only. 
These corrections have been applied based on a Monte 
Carlo simulation described below. Time intervals with 
background rates more than 4 times higher than the qui- 
escent rate were removed from the analysis to maximize 
the signal to noise in the RGS instrument. For the data 
selection cuts used here, the background comprises ~ 30% 
of the counts in the extracted spectrum. Approximately 
~ 70% of the 60 ks observation was used in this analysis. 

We detect two Fe L complexes dominated by Fe XXIV 
lines (10.6 A, 11.4 A) and O VIII Ly a (18.97A) above the 
bremsstrahlung continuum. The lines are broadened by a 
sharply peaked intensity distribution of ~ 1 arcminute 
FWHM. There are some non-statistical fluctuations, but 
all of these are much narrower than true emission lines 
from the cooling-flow would appear, and none of them are 


coincident with known strong astrophysical emission lines. 
The presence of the Fe XXIV ion suggests cool gas with a 
temperature between 1 and 3 keV, whereas O VIII could 
have contributions from both cool (~ 1 keV) and hot gas 
(~ 10 keV). The absence of other Fe L shell ions produces 
strong constraints on cooling-flow models as explained be- 
low. 

To fully address spectral-spatial coupling in the RGS, 
we compare the data directly to models by folding the 
spectral-spatial models through an instrument response 
Monte Carlo. The Monte Carlo takes wavelengths and sky 
positions of photons, and predicts cross-dispersion coor- 
dinates, dispersion coordinates, and CCD pulse-heights. 
The Monte Carlo includes characterizations of all known 
physical properties of the RGS instrument at the time 
of the writing of this paper, with the exception of a not 
yet calibrated O K instrumental edge of approximately 
20% at 23.5 A, which has moved to 18.7 A in our correc- 
tion of the redshift of Abell 1835 in Fig. 4. The analysis 
is then unaffected by the details of data selection cuts 
provided the same cuts are applied to the data and the 
Monte Carlo simulated photons. The cluster emission is 
distributed across the entire CCD array so there is no 
source-free region to estimate the background contribu- 
tion. Instead, we use a model based on the Lockman Hole 
(Revolution 70 and 73, lOOks observation) which charac- 
terizes the effects of low energy protons, the calibration 
sources, and the CCD detector noise. The background is 
smooth and is relatively flat across the wavelength band. 

In order to characterize spatial structure for use in the 
RGS Monte Carlo, we use a previous ROSAT observations 
for the density parameterization, which is consistent with 
the observed EPIC spatial distribution. We use a 0 model 
(King 1966, Cavaliere & Fusco-Fermiano 1976) to describe 
the density structure. The model has a core radius of 18 
arcseconds, 0 of 0.72. This characterization reproduces the 
RGS cross-dispersion event distribution. The exact spa- 
tial distribution is unimportant since it only affects the 
shape of the spectral line profiles and only a rough spa- 
tial characterization is needed. For the spectral models we 
use the coronal plasma MEKAL model of Mewe, Kaastra, 
and Liedahl (1995), the cold absorption model of Morri- 
son and McCammon (1983), and a spectral-spatial renor- 
malization of the MKCFLOW model for the cooling-flow 
emission measure included in the XSPEC vll.O package 
(Arnaud 1996). 

We select the spectral model for the Monte Carlo sim- 
ulation as a simple isothermal model folded through the 
instrument response based on the 8.2 keV measured tem- 
perature from the EPIC-MOS fit for the ambient gas. We 
also add a background component. We set all elemental 
abundances other than that of oxygen to 1/3 of solar val- 
ues and we set the iron abundance to 0.35 of solar values. 
The spatial distribution is given by the above mentioned 
beta model. We see that the agreement with the mea- 
sured continuum is quite good (see Fig. 4). The oxygen 
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Table 3. Cooling Flow Properties: RGS Monte Carlo* 


Parameter 

Isothermal 

Cooling- 

Cooling- 



Flow 1 

Flow 2 

^ambient (keV) 

8.2 

8.2 

8.2 

^Fe „„ 

0.35 

0.35 

0.35 

N H (10 2O cm -2 ) 

2.32 

2.32 

2.32 

0c(arcseconds) 

18 

18 

18 

0 

0.72 

0.72 

0.72 

r cool( k P c ) 


370 

370 


0.31 ±0.07 

0.19 ±0.04 

0.28 ± 0.07 

M(Af 0 yr- 1 ) 

*^nin( keV ) 

7V{J(10 2O cm" 2 ) 


1200 

1800 
2.7 ±0.2 


o +7 -° 
U '- 0.0 

q + 6.0 
U - 0.0 

Reduced \ 2 b 

1.08 

2.34 

0.92 

Reduced \ 2 c 

1.59 

3.46 

1.11 


* Parameters with no errors were held fixed 
b For the entire wavelength band (288 degrees of freedom) 
c For the 10 to 18 A region (100 degrees of freedom) 


line flux is well-matched by adjusting the abundance to 
0.31 ±0.07 of the solar value. The Fe XXIV complexes are 
above the expectation for the continuum and are directly 
attributable to the cooling- flow. 

Next, we attempt a multi-phase cooling flow model 
where the ambient gas is fixed at a temperature of 8.2 
keV and a fraction of the cool gas is expected to exist 
at all radii at several temperatures. Following the model 
of Johnstone et al. (1992), we model the cool component 
as having the same spatial distribution as the ambient 
component but only existing within the cooling radius. 
We have insufficient counts in the spectral lines to deter- 
mine the spatial distribution of the cool gas relative to the 
hot ambient gas. The emissivity of the cooling-flow gas is 
calculated from a conventional isobaric cooling-flow pre- 
scription, in which each temperature interval is weighted 
with the time the gas spends cooling through that interval. 
The resulting differential emission measure is then propor- 
tional to the inverse of the total radiative cooling function. 
We also allow for additional cold absorption applying it 
only to the cooling flow component at the redshift of the 
source, where the flux is reduced by (l*-e“ T )/rto account 
for self- absorption. This would apply if large deposits of 
cold byproducts of the cooling flow exist uniformly in the 
cooling-flow region. We fix the foreground column den- 
sity to the galactic value and apply it to both spectral- 
spatial components. The relative normalization between 
the two components is determined by the mass deposition 
rate (M) and the overall normalization is determined by 
central electron number density. 

To match the observed equivalent widths of the Fe 
XXIV lines (10.66 A, 10.62 A, 11.03 A, 11.18 A, 11.43 
A) , we require a mass loss rate of 1200 Mq yr -1 , but the 
model predicts significant line emission from several other 


iron ions that is not observed. This is seen quite easily 
from the blue model in the insert in Fig 4. In particu- 
lar, Fe XVII (15.01 A, 17.05 A, 17.10 A, 16.78 A) lines 
are most notably absent. Other ions such as O VII (21.6, 
21.8, 22.1 A) are not observed, but are not expected to be 
significant since the gas should cool very rapidly at lower 
temperatures. 

Previous ASCA studies found an emission line deficit 
at low energies (Fabian 1994), but this was modeled as uni- 
formly distributed additional absorption of ~ 10 21 cm" 2 
(Allen et al. 2000). Here we find no evidence for additional 
absorption uniformly distributed (TVjj < 6 x 10 2O cm 2 ). 
A larger uniform column density would have produced 
detectable absorption of the continuum seen through the 
cooling flow. We do not see this absorption so we model 
the spectrum by cutting off the temperature distribution. 
The EPIC fits to cooling flow models required either ad- 
ditional cold absorption of ~ 10 21 cm" 2 or a cut-off in the 
temperature emission measure distribution at ~ 2 keV. 

To characterize the spectrum empirically, we use the 
same cooling-flow model but put a lower temperature cut- 
off on the emission measure distribution of 2.7 keV as 
shown by the green model in Fig. 4. This cut-off is strongly 
constrained by the lack of Fe XXIII emission. The Fe 
XXIV 10.6 A complex is not well-matched because the 
true emission measure distribution probably falls more 
rapidly than implied by the isobaric cooling model. Fe 

XXIII contributes partly to the 11.4 A complex (Savin 
1996). 

If an additional cooling-flow component is added with 
no cutoff on the emission measure distribution, then its 
mass loss rate is less than 100 M 0 yr" 1 (90% confidence). 
The strongest constraint comes from the Fe XXIII, Fe 
XII, Ne X complex at 12.2 A. A smaller mass flow rate 
may be more consistent with the detection of molecular 
gas from CO emission, an indicator of star formation, of 
1.4 x 10 11 Mq. (A. Edge, private communication) Simi- 
larly Allen (1996) and Crawford et al. (1995) find optical 
spectra consistent with 100 M© yr" 1 . Another interpreta- 
tion of the data would be that a very small fraction (5%) 
of the cooling region is actually cooling and gas at in- 
termediate temperatures coexists, which accounts the Fe 

XXIV lines. 

4. Discussion 

The emission measure of the cooling-flow spectral model 
has to be artificially cut off at about 2.7 keV. This would 
be possible if the cooling time was reasonably long com- 
pared with the expected age of the cluster. There are three 
possibilities for this lack of emission. 1. It possible that 
there is an additional transport mechanism that has tem- 
pered the cooling and the becomes a frustrated cooling- 
flow quickly reheating the cool gas back to the the hot 
phase. 2. The gas has rapidly cooled below 2.7 keV and 
we therefore do not see it. 3. The cool gas may exist, but 
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Fig, 4. RGS spectrum of Abell 1835 and three models folded through the instrument Monte Carlo. The same corrections and 
data selection cuts have been applied to the data and the simulated photons. The red model is an isothermal 8.2 keV model. 
The blue model has an hot ambient 8.2 keV component and an isobaric cooling flow component. The green model is the same 
as the blue model but does not have emission below 2.7 keV. The details of each model are described in the text. The spectrum 
is corrected for redshift, exposure, and effective area. 


its emission is absorbed by cold material at the center of 
the flow. It is also important to note that the other XMM- 
Newton observations of clusters show very similar results 
suggesting that this is a generic problem in cooling-flow 
clusters (see Kaastra et al. 2000, Tamura et al. 2000). Fol- 
lowing, we discuss some possibilites. 

Underestimate of the Cooling Function: The pre- 
sented results depend somewhat on the cooling function 
assumed below ~ 1 keV, but the function would have to 
be underestimated by a factor of at least 10 to explain the 
results. However, for a gas at AjT e = l keV, we detect pri- 
mary radiative cooling directly in the X-ray emission lines 
in our band. From deep observations of coronal sources 
(Brinkman et al. 2000), we know that the emission spec- 
trum is accurately accounted for by the current spectral 
models, and hence the cooling function does not have a 
large error. Our results depend, however, on the metals 
IGM being approximately homogenously distributed. 

Resonant Scattering + Cold Absorption: This 
would be expected to be important in spectral lines with 
high oscillator strength where the optical depth, r ~ 10 
(e.g. 15 A Fe XVII could be r ~ 100). Resonance line pho- 
tons could then scatter severed times within the medium, 
which would appear to enhance the probability of photo- 
electric absorption. However, because resonance scattering 
of a single transition also has the effect of scattering the 
photon in frequency, the photoelectric absorption proba- 
bility is not enhanced by a large factor. Empirically, lines 


with smaller oscillator strength, which axe not subject to 
resonance scattering at all, are not observed (e.g. 17.1 A). 
The 10.6 A complex from Fe XXIV contains two medium 
oscillator strength lines which are still observed. 

Differential Cold Absorption: A column of 5 x 
10 21 cm~ 2 being selectively applied to the gas below 3 
keV would explain the observed spectrum. The limit of 
6 x 10 20 cm -2 on the column density being applied to the 
hot gas, however, would make this a very specific partial 
covering model. This could be reasonable, since the coolest 
gas occupies the smallest volume and is likely to be near 
the coldest material. The data, however, are inconsistent 
with the entire cooling volume being occupied by cold ma- 
terial. Further studies may help to constrain the mass and 
location of the cold material as well as the volume fraction 
of the cooling-flow. 

Non-Equilibrium Ionization: This effect would 
cause most ions to radiate at a lower temperature than is 
expected from collisional equilibrium models. (Edgar and 
Chevalier 1986). The recombination time scale for Fe XVII 
at 10 7 K is (cMie)” 1 « (10“ 12 cm 3 s” 1 x 10 -2 cm" 3 )' 1 « 
2 x 10 6 yr (Arnaud 1992). The cooling time, however, is 
« 10 9 yr, so this is unlikely to be important. 

Evolutionary Effects: The cooling time argument 
assumes that evolutionary effects are unimportant. A re- 
cent cluster merger, for example, would cause an overes- 
timate of the cooling radius. This seems unlikely since we 
would have to be fortunate enough to see the cooling-flow 
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just beginning. Another possibility is that there is secular 
mass injection throughout the lifetime of the cluster. This 
also seems unlikely since the mass of the IGM is larger 
than the mass of the galaxies themselves. 

Non-Isobaric Cooling-Flow: The cooling time es- 
timate is based on the assumption that the gravitational 
potential can be neglected. The cooling-flow will be non- 
isobaric if it is flowing into a region of a greater gravi- 
tational potential. This could apply if the gravitational 
potential does significant work on the cool gas (see e.g. 
Nulsen et al. 1982, Wise 1993). This would cause an over- 
estimate of the mass deposition rate and adds energy to 
the gas if the cooling-flow actually flows. 

Electron Thermal Conduction: Theoretical calcu- 
lations have disagreed about how much conduction is su- 
pressed by twisted magnetic fields (see e.g. Tao (1995) and 
Chandran and Cowley (1998)) and observationally it ap- 
pears to be suppressed (Ettori and Fabian 2000). Electron 
thermal conduction would have to be suppressed enough 
to allow for the observed cooling and temperature gradi- 
ents, but would have to play a role in limiting cooling. 
Note that cooling works best to balance cooling at higher 
temperatures. 

Magnetic Reconnection: If large scale magnetic 
fields were present in clusters, they might have the ef- 
fect of isolating cooling flow regions as well as reconnect- 
ing cooling-flow regions with the ambient gas (Soker and 
Sarazin 1990, Norman and Meiksin 1996). This would have 
the effect of recycling the cool material and delaying the 
cooling time. 

Turbulent Mixing: Begelman and Fabian (1990) 
suggested that turbulent mixing could result in an inter- 
mediate temperature layer between the hot gas and the 
cold condensed cooling-flow byproducts. If this process 
were highly efficient the cold gas could get recycled into 
an intermediate phase and the hot gas could be cooled 
rapidly by turbulent mixing within the cold layer. 

Dust Sputtering: Abell 1835 is an IRAS source 

(Allen 2000). This could imply that the dust is mixing 
with the cooling-flow and rapidly cooling the gas. 

The ideas above are only a sample of the possible mech- 
anisms for the absence of the 1 keV cooling-flow gas and 
more than one could be important. The observational evi- 
dence places strong constraints on cooling flow models and 
proper simulations of these effects are beyond the scope 
of this paper. With future XMM-Newton observations us- 
ing spatially-resolved spectroscopy and or further limits 
or detections of Fe XVII of many clusters, we expect to 
be able to distinguish between the alternatives. 

5. Conclusion 

We describe the spectral-spatial results of the observa- 
tion of Abell 1835 using the European Photon Imaging 
Cameras and the Reflection Grating Spectrometers on 
the XMM-Newton observatory. The observations are con- 


sistent with large amounts of cool gas surrounded by a 
kT e =8.2 keV component outer envelope. We find no evi- 
dence for excess uniform cold absorption above the galac- 
tic value of 2.32 x 10 2O cm -2 . The emission measure of 
cool gas below kT e =2.7 keV is much lower than expected 
in simple isobaric cooling flow models, and it indicates ei- 
ther differential cold absorption applied to the cool gas or 
that a more complicated description of the cooling process 
is required. 
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